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1. Introduction 

Numerical simulations based on lattice QCD allow us to calculate physical quantities in high 
accuracy []l|]. However, electromagnetic and isospin- violating effects are usually ignored. Taking 
into account these effects is important for more realistic predictions. In practice, it is pointed out 
that for light quark masses systematic errors from electromagnetic effects can be comparable with 
statistical one [Q]. 

An attempt was made to include electromagnetic effects in lattice QCD [H]. Pseudoscalar 
meson masses are calculated on the background of gluon and photon fields. n + — n° splitting is 
found to be 4.9(3) MeV at j8 = 5.7 with the Wilson action. Their result shows a good agreement 
with experiments. But their simulations were performed at a single lattice spacing with unimproved 
actions. The scaling violation of their result is not expected to be small. Another point is that the 
finite size corrections were estimated only with a model. It is desirable to evaluate finite size effects 
from a first principle calculation directly. 

In this work, we study electromagnetic mass splittings of mesons and extract quark masses 
in the quenched theory. 1 We employ an improved action combination, the renormalization 
group(RG) improved gauge action and the meanfield-improved clover quark action. Pseudoscalar 
and vector meson masses are computed at three lattice spacings and extrapolated to the continuum 
limit. Finite size effects are evaluated on lattices with the spacial size of L = 2.4 fm and L = 3.2 fm. 
Comparison of ra+ — rri\ as well as m+ — m° p with experimental values are made. Electromagnetic 
contribution to kaon mass splittings is also discussed. 

2. Method 

We generate SU (3) and U(l) fields and calculate quark propagators on the combined SU (3) x 
U(l) configurations. St/ (3) configurations are generated by the pseudoheat bath algorithm. For 
the U(l) gauge part we employ a non-compact Abelian gauge action. U (1) configurations are 
constructed of the Fourier transformed photon fields in the momentum space We fix the gauge 
to the Coulomb one in both SU (3) and £7(1) parts and use a gaussian smearing function. Meson 
masses are computed from quark propagators. In this calculation, we use only connected part 
of quark propagators. In the presence of electromagnetic fields, isospin is no longer a conserved 
quantity. Though the violation for m 2 K0 are expected to be small, disconnected diagrams may give 
measurable contribution. Estimate of the disconnected part is under investigation. Masses are 
obtained by x 2 fits to hadron correlators, taking account of correlations among different time slices. 
Statistical errors of hadron masses are estimated with the jack-knife procedure. 

Our simulations are performed at j3 = 2.187 (a p = 0.20 fm), j3 = 2.334 (a p = 0.16 fm) and 
j3 = 2.575 (a p = 0.11 fm) using 12 3 x 24, 16 3 x 32 and 24 3 x 48 lattices with the spatial extent L ~ 
2.5 fm. 16 3 x 24 lattices at j3 = 2.187 are also employed to examine the finite size effects. We take 
four hopping parameters corresponding to nips/my = 0.76 — 0.51. Configurations are generated 
independently at each mps/my. Measurements are carried out at each 100 sweeps. Our simulation 
parameters are summarized in Table [j]. These parameters are chosen so that they correspond to 
those of two-flavor full QCD data generated by CP-PACS. 

'For a similar work using Nf = 2 DBW2 gauge and domain wall fermion actions, we refer to Ref. [Q]. 
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Table 1: Simulation parameters 



j3 


2.187 


2.334 


2.575 


Size 


12 3 x 24(16 3 x 24) 


16 3 x 32 


24 3 x 48 


Nconf 


800(400) 


400 


100 



0.570 
0.568 
0.566 
0.564 
0.562 ■ 
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Figure 1: Volume dependence of charged pseudoscalar (left panel) and vector meson masses (right panel) at 
mps/mv — 0.62 on 12 3 x 24 lattice. For comparison, masses in pure QCD (open symbols) are also plotted. 



3. Simulation results 



We first check finite size effects in our results. Finite size corrections may be enhanced in the 
presence of electromagnetic fields because electromagnetic fields have a long interaction range. It 
is necessary to estimate magnitude of finite size corrections. In Fig. [I], we plot charged pseudoscalar 
and vector meson masses as a function of the spatial volume. The results obtained on 12 3 x 24 and 
16 3 x 24 lattices are mutually consistent within errors. We did not find any enhancement of finite 
size effects by electromagnetic fields. In the quenched approximation, L = 2.4 fm seems to be 
enough for meson calculations, even if there is an electromagnetic interaction. 

In order to extrapolate our results to the chiral limit, we fit a hadron mass as a function of 
quark masses and charges. We employ the following form for a chiral extrapolation of pseudoscalar 
meson masses. 

m PS = A Ps{q q + q-q) 2 + S PS {q q ,q q )(m q + m q ), (3.1) 
Bps{q q ,qq) = B PS +B PS (q q +q q ) 2 , (3.2) 

where Ap$,B PS ,B PS and K c in m q ,m q are fitting parameters. q q is an electric charge for a particle 
in units of e and q q is for an antiparticle. q = +2/3 for the up quark and q = —1/3 for the down 
quark are assigned. We use a quark mass defined through a vector Ward identity, 




(3.3) 
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Figure 2: Chiral extrapolations of pseudoscalar meson masses at /3 = 2.187. 



where K q is a hopping parameter for a quark and Kq is for an antiquark. fc^ are the corresponding 
critical hopping parameters. This extrapolation function, inspired by the chiral perturbation theory, 
only depends on the total charge and the quark mass combination. In the fits, correlations among 
several charge combinations are neglected for simplicity. We estimate the errors by the jackknife 



method. Strictly speaking, the quenched chiral logarithm term must be added to Eq. (p.ip. How- 
ever, the logarithmic curvature is not seen in our data within mps/my = 0.76 — 0.51. Smaller quark 
mass data are needed for a more precise extrapolation. We found Eq. ( fO| ) reproduce our lattice 
data well. As an example, the fit result at j3 =2.187 is presented in Fig. |2[ n + is heavier than 71°, 
as we see in nature. 

Vector meson masses are extrapolated as follows. 



m v = A v (q q , q-q) + B v (q q , q d ) (m q + niq) , 
Av(q q ,qq) =Al+A V 2 {q q + qq) 2 , 
Bv{qq,qq)=Bl+B v 2 (q q + qq) 2 . 



(3-4) 
(3.5) 
(3.6) 



Aq,A 2 ,Bq and B\ are fitting parameters. In contrast to the case of pseudoscalar mesons, mass dif- 
ference of charged and neutral vector mesons is found to be small. This tiny mass difference is con- 
sistent with the result calculated by the hidden local symmetry formulation, Am p = — 1 MeV 

We identify the physical point with experimental values of 71° and p° masses, M n o = 0. 1350 GeV 
and M p o = 0.7751 GeV. By solving Eq. Q3.1[ ) and Eq. ( ^4| ) using M^o^MpO, a sum of bare up and 
down quark masses at the physical point, (m u + m c j) p hy S , are determined. The lattice spacing a p is 
set by identifying M p o with my. Substituting (m u + md) p hys to Eq. ( p.l[ ) predicts n + mass. Sim- 
ilarly, Eq. (3.4) gives m p +. To determine up, down and strange quark masses independently, we 
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Figure 3: Lattice spacing dependence of quark masses with and without electromagnetic effects. 



also use experimental values of K + and Kq meson masses as inputs. (m u + m s ) p i iys is obtained from 
Mk + = 0.4937 GeV and (m^ + m s ) p h ys from Mk = 0.4976 GeV. Quark masses are renormalized 
using one-loop renormalization constants Z m and coefficients b m at /i = I /a [Q], 

m f^ = l/a)=Z m (l+b m T ^) r ^, (3.7) 

q \ UqJ Uq 

where uq is a tadpole factor, uq = (1 — 0.8412/jS) 1 / 4 . The MS quark masses at /I = I /a are evolved 
to /i = 2 GeV using the four-loop beta function. 

Continuum extrapolations are performed by linear fits to the data at three lattice spacings. 
Our preliminary results for light quark masses are m^ s {\i = 2 GeV) = 3.03(19) MeV, Tr^ s [\i = 
2 GeV) = 4.44(28) MeV, = 2 GeV) = 99.2(52) MeV. Fig. | illustrates lattice spacing de- 

pendence of quark masses with and without electromagnetic effects. We found electromagnetic 
contributions to the strange quark mass is 1%. In contrast to the case of quark masses, electro- 
magnetic mass splittings of pseudoscalar and vector mesons show mild lattice spacing dependence. 
Therefore, we employ constant fits for continuum extrapolations of electromagnetic mass split- 
tings. Our results are represented in Fig. f|. The obtained mass splittings are consistent with 
experimental values. In addition to n + — 7t° and p + — p mass difference, we can estimate elec- 



tromagnetic contribution to K + — K° mass difference by Eq. (3T). Using a constant continuum 
extrapolation again, the electromagnetic contribution to K + — K° mass difference is evaluated to be 
1.420(24) MeV, which is close to the value of Dashen's theorem 1.3 MeV [Q], rather than a model 
estimate 2.6 MeV But, our simulations are in the quenched approximation and the chiral log- 
arithm is neglected. Including dynamical quarks and chiral logarithm effects is needed for a more 
precise comparison with other calculations. 



4. Conclusions 

We calculated electromagnetic mass splittings using the RG-improved gauge action and the 
clover-improved Wilson quark action on the background of SU (3) x t/(l) fields in the quenched 
approximation. After chiral and continuum extrapolations, we found n + — n° and p + — p° mass 
differences are consistent with experimental values, and we can extract up, down and strange quark 
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Figure 4: 

value. 



Lattice spacing dependence of electromagnetic mass splittings. Star represents the experimental 



masses. We also confirmed that mass differences obtained with L = 2.4 fm are not shifted by finite 
size effects. An important future work is to include dynamical quark effects. Simulations using full 
QCD data generated by CP-PACS collaboration are ongoing. 
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